Giant magnetoresistance in ultra-small Graphene based devices 
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By computing spin-polarized electronic transport across a finite zigzag graphene ribbon bridging 
two metallic graphene electrodes, we demonstrate, as a proof of principle, that devices featuring 
100% magnetoresistance can be built entirely out of carbon. In the ground state a short zig-zag 
ribbon is an antiferromagnetic insulator which, when connecting two metallic electrodes, acts as a 
tunnel barrier that suppresses the conductance. Application of a magnetic field turns the ribbon 
ferromagnetic and conducting, increasing dramatically the current between electrodes. We predict 
large magnetoresistance in this system at liquid nitrogen temperature and 10 Tesla or at liquid 
helium temperature and 300 Gauss. 

PACS numbers: 



The term spintronics has been coined to refer to the 
interplay between spin polarization and electrical conduc- 
tance and is one of the major themes today in condensed 
matter and applied physics. A central concept in spin- 
tronics is that of giant magnetoresistance (GMR), discov- 
ered originally in layered structures alternating magnetic 
and non-magnetic transition metals The resistance of 
the whole structure undergoes a large increase when the 
relative orientation of the magnetization in adjacent lay- 
ers goes from parallel to antiparallel, provided that the 
non-magnetic layers are thinner than the spin relaxation 
length. This phenomenon represents the paradigm un- 
der which commercial devices such as magnetic reading 
heads operate nowadays. 

Here we propose a new type of magnetoresistive device 
which makes use of the remarkable electronic properties 
of zig-zag graphene ribbons [2, 3]. Its operational princi- 
ple is similar to that found in conventional GMR layers, 
with the difference that is entirely based on carbon. Our 
proposal is motivated by the spectacular progress in the 
fabrication of high mobility graphene based field effect 
transistors [3, H, 0] and by the recent developments in 
the fabrication of graphene ribbons 0, S, E3] with top- 
down techniques, as well as in the sinthesys of graphene 
ribbons 11 1 . 



The electronic structure of infinite graphene ribbons 
has been studied thoroughly. Idealized graphene rib- 
bons, with boundaries parallel to the crystallographic 
directions and boundary carbon atoms passivated with 
a single hydrogen atom, fall into two categories: arm- 
chair and zig-zag. Within the simplest one orbital tight- 
binding description armchair ribbons can be either 
semiconducting or metallic, depending on their width, 
whereas in the case of zigzag ribbons two flat bands, as- 
sociated with edge states, lie at the Fermi energy. These 
edge flat bands favor the appearance of magnetization on 
the edges when electron-electron repulsion is included in 
the calculation, either with a Hubbard model 3, [T2I E3] 
or with density functional theory (DFT) [3, EH, Ei[ ■ In 



the ground state the respective magnetization direction 
of the edges is antiparallel, and a gap opens in the 
band structure 12, 3, EE EH- This is the ground state. 
Slightly above in energy, the parallel magnetic configu- 
ration is conducting. Application of either a magnetic 
field or a transverse electric field [15] can make the ferro- 
magnetic configuration more stable. Here we explore the 
former. 

Both for conceptual and practical reasons our pro- 
posal is based on finite length L graphene ribbons. For 
L = 00, the long range order predicted by mean field cal- 
culations is not robust, due to the proliferation of spin 
wave excitations of energy L~ n , with n = 1 or n = 2 
for both antiferromagnetic and ferromagnetic alignments, 
respectively fl7j . In short ribbons, in contrast, there is a 
gap for spin waves. Quantum fluctuations between the 
manifold of equivalent mean field ground states would 
not alter the conducting properties of the system. Addi- 
tionally, recent DFT calculations [l|| show that monohy- 
drogenated zigzag ribbons are not the most stable edge 
configuration unless the hydrogen density is very small. 
Whereas this poses a severe problem for the chemical sta- 
bility of infinite monohydrogenated zigzag ribbons, short 




FIG. 1: Atomic structure of the zigzag ribbon with length 



12 and width N v 



6 attached to semi-infinite elec- 



trodes. The unit cell of a zigzag ribbon is highlighted. 
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magnetic ribbons, as well as other magnetic polycyclic 
aromatic hydrocarbons [l9j| . might be more stable. 

Here we show that the spin driven metal insulator tran- 
sition, predicted so far for infinite ribbons, is still present 
in short ribbons attached to conducting electrodes. Using 
a well-established methodology[2(| extended to account 
for electron-electron interactions in a Hubbard model, 
we study both the magnetic and transport properties of 
a system in which two conducting graphene electrodes 
are coupled through a finite length graphene ribbon. 

Statement of the problem and methodology.- We study 
the system depicted in Fig. 1: A short zig-zag graphene 
ribbon attached to two conducting graphene electrodes. 
In the case pictured in the figure we have chosen the elec- 
trodes to be semi-infinite metallic armchair graphene rib- 
bons of width W . Both the device and the electrodes are 
described with a Hubbard model for the 7r z orbitals with 
first-neighbor hopping t and on-site repulsion U, solved 
at a collinear mean field level0, [H [H HI], El O This 
approach is known to capture the main features of the 
ab-initio calculations, both for finite 21 1 and infinite [l2| 
graphene systems. Thus, the mean-field Hamiltonian 
reads 



I, I'. a 



where c\ a creates an electron at the -k z orbital of atom 
I with spin a, ni y(T = c\ a ci a is the occupation operator. 
Since the mean fields (n-/,o-) depend on the eigenstates of 
the mean field Hamiltonian, this defines a self-consistent 
problem which is solved by numerical interation. The 
converged solutions define a one body Hamiltonian for 
the electrons in the structure. 

In this work we have to solve the self-consistent prob- 
lem for an infinite system without translational invari- 
ance. This is done using the partition method|20|. We 
split the system in three sectors, the left and right elec- 
trodes and the central region. The Hamiltonian of Eq. 
(HJ) reads: 



H = H L +H R + Hc + V LC + V lr 



(2) 



where 1~1{l,r,c) are the mean field Hubbard Hamiltoni- 
ans of the left and right electrodes and the central region 
and V(ij j L)c describes the hopping between the central 
region and the electrodes. Sufficiently away from the 
central region, the mean field Hamiltonian of the elec- 
trodes is identical to that of an infinite ribbon. In the 
case of the armchair electrodes, the effect of the Hub- 
bard interactions in the charge neutrality point is a rigid 
shift of the bands without the appearance of magnetic 
moment and keeping their metallic character. The first 
step is the determination of the surface Green function of 
the semiinfinite electrodes, g^L, R ) [20]. This requires the 
solution of a self-consistent Dyson equation. Once this is 



done, the Green function of the central region reads: 

G c (E) = [EI-Hc-^l(E)^J: r (E)}- 1 (3) 

with Z {L , R) (E) = V (L , R) g (LM) (E)vl L R) 

This expression is a functional of the expectation val- 
ues (ni M ), through the central Hamiltonian Tic- The 
Green function yields the density of states projected over 
the orbital ~k z with spin a sitting in the atom I in the 
device, is given by: 

p(E,I,a) = --ImTr((Ia\G c (E)\Ia)) (4) 

Here (Ia\Gc{E)\Ia) is a matrix of size 2N, with N being 
the number of atoms of the central region. In turn, the 
expectation values of the spin density arc given by 



(ni,*) = 



p{E,I,a)dE 



(5) 



The magnetic moment in a given atom I is defined as 

UlJ = <"'.T>-<"*a) 

Equations ((3|) and |5]) define a self-consistent problem 
which is solved by numerical iteration. The solution pro- 
vides a mean field description of the central region at- 
tached to the electrodes. Within this formalism the Lan- 
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FIG. 2: (a,b) Bands for infinite ribbon with N y = 12 with 
AF (a) and FM (b) ground states. (c,d) Corresponding DOS 
of the AF (c) and FM (d) infinite ribbons, (e) DOS of the 
AF finite ribbon, (f) DOS of the FM finite ribbon. 

dauer conductance is given by G — ^-TrT(Ep) with: 

T(E) = Tr (t l (E)G c T r {E)G ] c {E)) (6) 
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with Y L {E) = i {p {ttR) {E) - V\ LtR) {E)) 

Results.- The size of the zig-zag ribbons is defined by 
two integer numbers, N x , the number of unit cells of the 
ribbon, and N y the number of zigzag rows in the unit cell. 
Thus, the length of the ribbon is L = N x a, with a = 2.42 
A and the width of the ribbon is given by W = ^-N y a. 

We verify first that magnetic solutions appear in short 
zigzag ribbons connected to conducting electrodes. As 
in the case of both isolated finite ribbons [3] and infi- 
nite ribbons [15|. we obtain two kind of solutions, ferro- 
magnetic and antiferromagnetic, depending on the initial 
guess in the iterative self-consistent procedure. We see 
how only the edge atoms of the zigzag ribbon are mag- 
netic and how the size of the edge moments is larger away 
from the interface. In the case of the AF solution, the 
largest magnetic moment is mj = 0.13 for U = t = 2.7eV 
almost identical to the result of infinite ribbons. In Fig. 2 
we compare the electronic structure of the finite size con- 
nected zigzag ribbons with the case of infinite ribbons, 
both for ferromagnetic (FM) and antiferromagnetic (AF) 
configurations. The electronic structure of the infinite 
ribbons is calculated taking advantage of crystal invari- 
ance and making use of the Bloch theorem [l2j. Whereas 
AF infinite ribbons have a gap Ao 12, Id | and zero den- 
sity of states (DOS) at Ep (hereon set to zero), two bands 
cross Ep for the FM ribbons, resulting in a finite density 
of states at Ep. The same trend is observed in finite rib- 
bons. The DOS of the AF configurations, summing over 
the atoms of a unit cell (see Fig. 2), presents a pseudo- 
gap. The DOS becomes more similar to the one of the 
infinite case both for longer ribbons and for unit cells in 
the middle of the ribbon. The small DOS inside the gap 
arises from the coupling to the electrodes. In contrast 
to the AF solution, the DOS of the FM solutions, again 
summing over the atoms of a unit cell (see Fig. 2), is 
that of a conducting system. 
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FIG. 3: (a) Spin resolved transmission for the AF infinite 
ribbon with N x = 12 and iV y = 6. (b) Same for the FM 
solution. 



From the DOS we anticipate a strong dependence of 



the conductance on the magnetic configuration. Our cal- 
culations confirm this expectancy. In Fig. 3 we show the 
conductance curves for the N x = 12, N y = 6 ribbon. In 
the left panel we show how the transmission of the AF 
solution for E = Ep is highly reduced, compared to the 
FM solution, shown in the right panel. This is the central 
result of our work: The conductance is strongly depen- 
dent on the relative orientation of the magnetization of 
the edges of the zigzag ribbon bridging the electrodes. 
In the ground state the ribbon is antiferromagnetic and 
the conductance is small or vanishes for sufficiently large 
N x . In Fig. 4 we plot the conductance at Ep, both for 
AF and FM configurations, as a function of the ribbon 
length L = aN x for two widths. Only for very short rib- 
bons the conductance in the AF case is not vanishingly 
small, due to tunneling. The conductance of the FM 
state is always close to 2e 2 /h, value associated with the 
two bands, one per spin channel, that cross the Fermi 
energy. This conductance value is sensitive to the ac- 
tual form of the graphene electrodes, but the vanishingly 
small values for the AF state are not. 

As in standard GMR devices, the application of a mag- 
netic field can force the ferromagnetic solution to be the 
ground state, resulting in a dramatic increase of the con- 
ductance. We define the magnetoresistance (MR) of each 
device as 



MR: 



Raf — Rfm 
Raf + Rfm 



x 100 



G 



FM 



Gai 



Gpm + Gaf 



x 100, (7) 



where R = ^ is the resistance and G = j-T(Ep) is 
the conductance calculated with the Landauer formula. 
In Fig. 4 we plot the MR as a function of the ribbon 
length for ribbons of width N y = 6 and N y = 3. We 
see that, except for ribbons shorter than lnm, the MR 
saturates to 100%. The origin of the MR proposed here is 
different from that discussed by Brey and Fertig[l^] and 
that of Kim and Kim[26| which require ferromagnetic 
electrodes. Our proposal is more similar to the original 
GMR experiments with current in plane in exchanged 
coupled multilayers [l|. 

We now discuss the operational limits of our proposed 
device. There are three relevant magnitudes to be con- 
sidered. First, the energy difference between the AF and 
the FM states, which depends both on the width of the 
ribbon and the length: 



A(N x ,N y ) ~ N x 



30 



3/2 



(8) 



This expression has been obtained from the fitting 
to DFT calculations by Pisani et al.[l6| Second, the 
temperature-dependent spin correlation length £ (T) , over 
which magnetic order is lost [T3| • Finally, and neglecting 
the influence of the magnetic field on the orbital part of 
the wavefunctions, the critical switching magnetic field is 
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FIG. 4: (a) Conductance and (b) magnetoresistance (Eq. [7J) 
at the Fermi energy as a function of the ribbon length for two 
ribbon widths, both for FM and AF configurations. 



ribbon joining metallic graphene electrodes. The conduc- 
tion properties of this device change dramatically as the 
relative orientation of the magnetic edges of the ribbon 
go from parallel to antiparallel relative orientations upon 
application of a magnetic field. Even if from the oper- 
ational point of view there are still challenges regarding 
room-temperature performance in small magnetic fields, 
a proof of principle has been demonstrated. From a more 
fundamental point of view, low-temperature experiments 
showing a drastic change in the resistance on applying 
strong magnetic fields would unambiguously signal the 
existence of magnetism in zig-zag graphene ribbons. 

We acknowledge fruitful discussions with D. Soriano. 
This work has been financially supported by MEC- 
Spain (Grant Nos. MAT07-67845 and CONSOLIDER 
CSD2007-0010). 



defined through the equation 

gf i B M T B* = A (9) 

where Mr = J2i m i 1S the total spin density of the FM 
configuration, g — 2 for graphene and \xb = 0.058 meV 
T '. 

The temperature determines the minimal value A m i n 
below which thermal fluctuations will make the AF un- 
stable. It also determines the maximal value for the 
length of the ribbon L max above which magnetic order 
is lost. By choosing L — L max , we guarantee the max- 
imum possible width W max for the ribbon through Eq. 
|S) This, in turn, provides the minimum critical switching 
field B* through Eq. [5J For instance, at room tempera- 
ture, A m i n w 26 meV and £(300) = L max w 0.8nm. This 
yields W max « 1.6nm and B* « 200 T. For liquid nitro- 
gen temperatures (75 K) A m in ~ 6 meV, L max ~ 4nm, 
W m ax ~ 12 nm, and B* w 10 T. At He liquid tem- 
peratures (4 K) A min « 0.35 meV, L max ss 80nm, 
W max « 600 nm, and B* m 0.03 T. Using the extrap- 
olation formula of Son et aL[14| for the ribbon transport 
gaps a ribbon of W w 600 nm has a gap of 1.55 meV, 
which still larger than kT and standard bias voltages. 
In summary, there is a trade-off between the minimum 
temperature below which the AF ground state becomes 
unstable and the critical switching field. This trade-off 
is controlled by the spin correlation length. At high tem- 
peratures B* becomes prohibitively large. At liquid He 
temperatures, the critical field is bounded within reason- 
able ranges easily attainable in the lab. 

Also the width of the electrodes does not need to be 
equal to the width of the zigzag channel. This would only 
introduce some additional scattering in the FM state, but 
never compromising the high conductance values. 

In conclusion, we propose an ultra-small and chem- 
ically simple magnetoresistive device based on a zigzag 
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